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ABSTRACT 

Trojan asteroids are minor planets that share the orbit of a planet about the Sun 
and librate around the L4 or L5 Lagrangian points of stability. They are impor- 
tant solar-system fossils because they carry information on early Solar System 
formation, when collisions between bodies were more frequent. Discovery and 
study of terrestrial planet Trojans will help constrain models for the distribu- 
tion of bodies and interactions in the inner Solar System. Since the discovery 
of the first outer planet Trojan in 1906, several thousand Jupiter Trojans have 
been found. Of the terrestrial planets there are four known Mars Trojans, and 
one Earth Trojan has been recently discovered. We present a new model that 
constrains optimal search areas, and imaging cadences for narrow and wide field 
survey telescopes including the Gaia satellite for the most efficient use of tele- 
scope time to maximize the probability of detecting additional Earth Trojans. 

Key words: methods: numerical - methods: observational - minor planets, 
asteroids: general - planets and satellites: general - celestial mechanics 



1 INTRODUCTION 

Trojan asteroids are minor planets that share the orbit of 
a planet about the Sun and hbrate around the L4 and L5 
Lagrangian points of stability. The L4 and L5 points are 
60° ahead and behind, respectively, the planet in its or- 
bit. Trojans represent the solution to Lagrange's famous 
triangular problem and appear to be stable on long time- 
scales (100 Myr to 4.5 Gyr) (Pilat-Lohinger, Dvorak & 
Burger 1999; SchoU, Marzari & Tricarico 2005) in the N- 
body case of the Solar System. This raises the question 
whether the Trojans formed with the planets from the 
Solar nebula or were captured in the Lagrangian regions 
by gravitational effects. Studying the Trojans provides 
insight into the early evolution of the Solar System. 

Since the discovery of the first Trojan in 1906 
(Nicholson 1961) several thousand more have been found 
in the orbit of Jupiter. Of the terrestrial planets, four 
Trojans have been discovered in the orbit of Mars. There 
have been some attempts at searching for Earth Trojans 
(ET) (Dunbar & Helin 1983; Whiteley & Tholen 1998; 
Connors et al. 2000). Examination of a sky area of ap- 
proximately 0.35 deg^ by Whiteley & Tholen (1998) re- 
sulted in a crude upper limit on population estimates of 
~ 3 objects per square degree, down to i? = 22.8. The 
subsequent search by Connors et al. (2000) covered ~ 9 
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deg , down to i? = 22, with no ETs detected. Recent ex- 
amination of data from the WISE satellite has resulted 
in the discovery of the first known ET (Connors, Wiegert 
& VeiUet 2011). 

The earlier non-discovery may be attributed to a 
lack of observations targeting regions in which ETs could 
be found. For example, the programmes to search for 
Near-Earth Asteroids (NEA) generally employ observa- 
tion strategies which survey regions having a higher prob- 
ability of containing potential Earth impactors. Typically 
these regions are near the plane of the ecliptic and do 
not specifically target the entire region of stability for 
ETs. The Earth- Sun-ET geometry requires observations 
at low altitudes with limited observing time, restricting 
the ability to make multiple or follow-up observations. 
This geometry also results in low apparent motion rela- 
tive to field stars as the orbit arc is at an oblique angle 
to an observer on Earth. Other reasons for non-detection 
include small population and orbital inclinations outside 
the plane of the ecliptic. Asteroids which could be Tro- 
jan candidates would not be flagged for further study 
by routine surveys because their apparent motions would 
not match the parameters of the survey for follow-up. 

Models (Mikkola & Innanen 1990; Tabachnik & 
Evans 2000a) indicate ET stable orbits have some incli- 
nation to the plane of the ecliptic, in the range 10° to 45° . 
Subsequent simulations (Morais & Morbidelli 2002) pre- 
dict 0.65 ±.012 ETs with diameter > 1 km and 16.3 ±3.0 
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Table 1. Absolute and apparent magnitudes at the classical 
L4 and L5 Lagrangian points. 



Type 


Albedo 


Diameter 


Abs. mag. (H) 


App. mag. (V) 


S-type 


0.20 


1.0 km 


17.37 


18.9 






100 m 


22.36 


23.9 


C-type 


0.057 


1.0 km 


18.73 


20.3 






100 m 


23.73 


25.3 




90 70 50 

Longitude (deg.) 



Magnitude 

119,5 
19.3 
19.1 
16.9 
18.7 
18.5 
18.3 
18.1 
17.9 



asteroids with diameter > 100 m. The existence of the 
co-orbital asteroid 3753 Cruithne (Wiegert, Innanen & 
Mikkola 1997), and the recent discoveries of the horse- 
shoe orbiter 2010 SOie (Christou & Asher 2011) and 
L4 Trojan 2010 TK7 (Connors, Wiegert & Veillet 2011), 
confirm the possibility of objects orbiting the Sun in 1:1 
mean motion resonance with Earth. The question is how 
many, and what size are the bodies that share Earth's 
orbit as Trojan asteroids. 

Assuming Trojans are approximately evenly dis- 
tributed between the L4 and L5 regions then just one 
of these fields could be examined with the expectation 
that the other region would yield a similar result. Ifow- 
ever it should be considered that, given the short window 
of opportunity for making ground-based optical obser- 
vations near morning (L4 region) and evening twilight 
(L5 region), both regions should be surveyed. With this 
in mind, it is important to find the optimal strategy to 
maximize sky coverage and probability of detection given 
those time constraints. 

This paper employs a model probability distribution 
which we use to constrain optimal search areas and imag- 
ing cadences for efficient use of telescope time while max- 
imizing the probability of detecting ETs. 



2 MODEL 

Existing models for ETs provide estimates on popula- 
tions but not composition. This is an important factor 
considering that albedo depends on composition. In the 
most likely case, ETs may be of S-type (silicaceous) sim- 
ilar to NEAs and inner Main Belt asteroids, or they may 
be of C-type (carbonaceous). This influences the detec- 
tion limit as S-types have an albedo p„ = 0.20 and C- 
types have an albedo pv = 0.057 (Warner, Harris & 
Pravec 2009). Calculations of absolute and apparent mag- 
nitudes for the classical L4 and L5 Lagrangian points, 
from Tedesco, Cellino & Zappala (2005) and Morals & 
Morbidelli (2002), are shown in Table 1, neglecting at- 
mospheric extinction. Assuming S-types as the dominant 
class in the inner Solar System the apparent magnitude 
for an ET of 1 km diameter varies between V = 17.9 to 
V = 19.5 (Fig. 1), depending on the geometry across the 
field (Fig. 2). An ET of 100 m diameter varies in mag- 
nitude between V — 22.9 to = 24.5 across the field in 
the same fashion. 

A synthesis of the stable orbit inclination model 
(Morais & Morbidelli 2002) and heliocentric longitude 
model (Tabachnik & Evans 2000b) was used to identify 
probability regions where bodies are most likely to be 



Figure 1. Apparent magnitude of a 1 km Earth Trojan by 
Heliocentric Longitude and Latitude. The figure shows varia- 
tion from V = 17.9 at the nearest part of the field (to Earth) 
to V = 19.5 at the farthest part. 




Figure 2. Perspective illustration of Earth Trojan (L4) tar- 
get field. The field ranges from Heliocentric longitude (A) 30° 
to 130° and latitude (/3) 10° to 45°. A complementary field 
exists in the trailing Lagrangian L5 region. This illustration 
represents the field through which a body will pass during its 
orbit. 



(Fig. 3). The ET fields (Fig. 2) bounded by upper and 
lower inclination limits (FWHM) of 10° < /3 < 45° in- 
clude ~ 74 per cent of bodies (Morais & MorbideUi 2002). 
The heliocentric longitude limits (FWHM) of 30° < A < 
130° (L4 region) and 240° < A < 340° (L5 region) in- 
clude ~ 45 per cent and ~ 40 per cent of bodies respec- 
tively (Tabachnik & Evans 2000b). The regions bounded 
by these limits enclose ~ 63 per cent of projected bod- 
ies. The recently-discovered ET, 2010 TK7 (Connors, 
Wiegert & Veillet 2011), has inclination i = 20.88° and 
mean longitude within the limits for the L4 region. 

The heliocentric solid angle of each region is 0.931 sr 
(3056 deg^). Calculation of the geocentric solid angle, 
necessary for Earth-based observations, requires a trans- 
formation from the heliocentric reference. A numerical 
integration is performed using the standard solid angle 
integral JJ^ (r • n) /r^dS (where r is the distance vector 
from the geocentre to the surface dS, n is the unit vec- 
tor normal to that surface, and r — |r|) to determine 
the solid angle subtended at points other than the cen- 
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Figure 3. Normalised probability contour for Earth Trojan 
bodies by Inclination and Heliocentric Longitude (degrees). 
The figure shows peak detection probabilities for longitudes 
consistent with the classical Lagrangian points but that bod- 
ies, while co-orbital with Earth, are unlikely to be co-planar. 

tre. This integration enables calculation of sky area^ for 
the heliocentric surface, from the geocentre for an Earth- 
based observer or for a space-based instrument such as 
the Gaia satellite which will be located near the Earth's 
L2 Lagrangian point (Mignard et al. 2007). The calcu- 
lated geocentric solid angles are 0.3948 sr (1296 deg^) 
and 0.4451 sr (1461 deg^) for the L4 and L5 regions re- 
spectively. From the L2 Lagrangian point, for the Gaia 
satellite, the solid angles are 0.3855 sr (1265 deg^) and 
0.4310 sr (1415 deg^). The variation between regions is 
due to the slight differences in positions relative to Earth. 

2.1 Telescope surveys 

Having identified the ET fields we consider the ability 
of current and proposed wide-field telescopes to survey 
these regions. It is clear from the size of the regions 
(>1000 deg^) that a widefield survey telescope would be 
needed for any attempt at surveying the entire region. 
However even existing telescopes with the widest fields 
of view would barely be able to survey the entire region 
once in a single day due to observing constraints imposed 
by the Earth-Sun-ET geometry. 

We also consider the implications for a space-based 
telescope with the impending launch of the Gaia satel- 
lite. Gaia will not be subject to the same constraints as 
ground-based telescopes. Although Gaia has a relatively 
small FOV it will operate in a continuous scanning mode. 
It will image all of the sky down to a Solar elongation of 
45° (Mignard et al. 2007), but without the observational 
limitation of local horizon and airmass. However, this ad- 
vantage is mitigated by its limiting magnitude of 1^ = 20. 

The continuous-scanning operation mode of Gaia 
leads to the concept of a strategy of observing a swath 
of the region, a range in Right Ascension and Declina- 
tion defining a particular sub-region. Observations can 
be made of that defined region, repeated twice per week, 

^ The Python code used to calculate these solid angles is avail- 
able from the author on request. 




Figure 4. Observing a defined range in Right Ascension and 
Declination, with Earth's revolution about the Sun, results in 
the entire field being imaged as it passes across the observed 
region. 

and making use of Earth's revolution about the Sun to 
progressively survey the target field. A 10°-wide swath 
(Figure 4) would be imaged in minutes by a survey tele- 
scope, as shown in Table 2. The region is redefined at 
monthly intervals to recommence the survey. Thus each 
month the target field passes through the observed re- 
gion. In this fashion the entire field is resampled over a 
period of one year. This can be done at the beginning 
and end of the night, requiring minimal telescope time, 
before and after the primary science missions. 

Observations of ETs are time-limited by the geom- 
etry of the Earth-Sun-ET positions. There exist specific 
constraints particular to the geographic location of a tele- 
scope. For example Northern Hemisphere telescopes will 
have access to the entire Northern part of the ET region 
described in Figure 2 while Southern Hemisphere tele- 
scopes will have access to only that portion visible above 
the local horizon. Some seasonal variation will occur with 
the change in relative orientations throughout the year, 
as noted in Whiteley & Tholen (1998). The amount of 
variation will depend on geographic location. 

By the end of evening twilight when it becomes pos- 
sible to survey the L5 region, assuming (for convenience) 
the Sun's altitude is —20°, the classical Lagrangian point 
has an altitude of about 40°. Hence atmospheric extinc- 
tion on detection limits must be considered. Because 
there are many contributing factors to atmospheric ex- 
tinction, reference to the "average" approximation given 
by Green (1992) suggests the extinction at 40° altitude 
is 0.44 mag and at 25° altitude is 0.66 mag for a site at 
sea level. 

The stable orbits of the ETs have been modelled to 
have an inclination range of 10° - 45° (Morals & Mor- 
bidelh 2002). The heliocentric longitude ranges (FWHM) 
are 30° - 130° and 240° - 340° for L4 and L5 respectively 
(Tabachnik & Evans 2000b). Since the modelled orbits 
are inclined at some angle to the ecliptic the survey could 
be carried out in the upper or lower ecliptic latitudes (/3) 
at 10° ^ /3 s; 45° and -10° < ;fl < -45°. By surveying 
fields in those latitudes the transverse motion per unit 
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Table 2. Comparison of different survey telescopes showing the required time to survey the entire field compared to an optimal strategy 
of observing a 10° swath and using Earth's revolution to sweep across the field. 



Telescope 


Limiting mag. 


Exposure 


FOV 


Entire field 


10° 


swath 


Instrument capabilities 










(FOVs) 


(Time) 


(FOVs) 


(Time) 




Catalina 


y ~20 


30 s 


8.0 deg2 


162/183 


1.35h/1.5h 


12-18 


6-9 min 


(Drake et al. 2003) 


PTF 1.2 m 


R ~ 20.6 


60 s 


8.1 deg2 


160/180 


2.7h/3h 


12-18 


12-18 min 


(Law et al. 2009) 


Pan-STARRS 


R ~ 24 


30 s 


7.0 dcg2 


185/209 


1.5h/1.7h 


13-20 


7-10 min 


(Jcdickc ct al. 2007) 


LSST 


r ~ 21.7 


30 s 


9.6 dcg2 


135/152 


l.lh/1.3h 


10-15 


5-8 inin 


(.loiios et al. 2009) 


Gaid 


I" - 20 


:-!!). (i si 


().()!)- 


():?/7i 


:?781i/ililit 






(Mignard c-l al. 2007) 



^Gaia will operate in a continuous scanning mode where the CCD array will be read out at a rate corresponding to the angular rotation 
rate of the satellite (6h period). The FOV value represents the number of rotations by Gaia. Gaia's specific precession parameters are not 
considered so values should be considered as representative. 

JApproximate time to complete sufficient rotations to scan across the entire field. The actual time spent scanning within the field will be 
a fraction of this value. 



time is reduced. To an Earth-based observer an ET will 
appear to oscillate in an up-and-down motion about the 
plane of the ecliptic during its orbit about the Sun. The 
exhibited transverse motion will appear sinusoidal, with 
the apparent transverse motion reaching minimum at the 
extremes of its heliocentric latitude and reaching maxi- 
mum when crossing the ecliptic plane. 

It is impossible, a priori, to predict the direction of 
motion for the ET. Consequently it is not possible to ap- 
ply a tracking offset for the apparent motion to increase 
the detection limit. A suitable alternative is to observe 
with the intention of imaging the ET at maximum he- 
liocentric latitude when the apparent motion is at mini- 
mum. As it will then appear to be a slow-moving object, 
follow-up becomes difficult in the same night. However, 
repeat imaging can be performed on subsequent days 
rather than trying to acquire follow-up images during the 
same session. This approach also has the advantage of in- 
creasing the available observing time for a target field, if 
necessary. 

While this delay between follow-up images intro- 
duces other variations from such things as changes in 

atmospheric conditions and seeing, this could be com- 
pensated for by image convolution. Some telescopes are 
implementing image processing systems designed specif- 
ically for asteroid detection (e.g. Pan-STARRS-|-MOPS 
- Moving Object Processing System) where the asteroid 
is an apparently stationary transient because it is a very 
distant slow moving object (Jedicke et al. 2007). This 
method can be applied to nearer objects which are ap- 
parently slow moving as a result of the particular Earth- 
Sun-MP positions at the time of observation. 



3 RESULTS 

Surveys of the entire field within the chosen limits are 
impractical on telescopes with small FOV but arc rea- 
sonable on survey telescopes with sufficient FOV to ac- 
complish the task in a single night, such as Catalina or the 
Large Synoptic Survey Telescope (LSST)^. In these cases 
it may be possible to survey the regions defined by the in- 
clination limits at intervals of approximately two months. 

The LSST is still in the development phase (www.lsst.org). 



This requires the fields to be surveyed twice within a few 
days. These surveys can be conducted at the beginning 
(L5 region) and end (L4 region) of the night. Although 
this is possible, it may be impractical as this occupies a 
significant amount of telescope time on those nights. 

As the field spans 90° in longitude the intervals could 
nominally be three months. This increases the risk of 
missing low-inclination ETs completely through not hav- 
ing yet entered the field or of having just exited the field. 
It also increases the risk of detecting ETs shortly before 
exiting the field towards the ecliptic, possibly requiring 
follow-up observations in a less favourable orientation to 
the local horizon. By observing at intervals of two months 
a small amount of oversampling is introduced and the risk 
of missing an ET is reduced. On the assumption that 
only the Northern or Southern field could be surveyed 
depending on the location of the telescope, the program 
length would be 1 year, the length of the orbital period 
of a Trojan as it would pass through this field once per 
revolution. 

It is possible to conduct a survey in the ecliptic plane 
between latitude and the lower limit of the field to 
search for Trojans as they cross the ecliptic. This reduces 
the necessary time per-session to survey the region as the 
ecliptic region is smaller than the ET field. The length of 
this program is reduced, requiring six months. Any Tro- 
jans detected in that six month period would be crossing 
either to the North or to the South. However the time 
saving per session is only about 30 per cent and requires 
more frequent sampling since, in this region, the ET will 
have a higher apparent motion. This could outweigh any 
benefit in the reduced per-session length as the total time 
requirement is not significantly different to that required 
to survey the entire field over the period of one year. 

Attempting to survey the entire field in a single ses- 
sion is challenging. The desirability of minimising the 
time per session is our key objective. Observing a swath 
of sky in a particular range of Right Ascension and Dec- 
lination, and using Earth's revolution about the Sun to 
survey the field £is it crosses the observed region, achieves 
this aim. This approach requires minimal time each ses- 
sion, i.e. two sessions per week during the course of one 
year. The observed region is redefined at monthly inter- 
vals to recommence sampling of the field. Some oversam- 
pling will occur as the field crosses the observed region. 
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This results in nights lost due to adverse weather condi- to examine the potential for discovery of these lEOs dur- 
tions becoming less critical to the overall program. ing an ET search. 



4 SUMMARY AND FUTURE WORK 

Despite the thousands of known Jupiter Trojans, of the 
terrestrial planets a mere handful of Trojans have been 
discovered. There are only four known Mars Trojans, and 
the first ET has only been very recently discovered (Con- 
nors, Wiegert & Veillet 2011). Simulations (Morais & 
Morbidelli 2002) have predicted the existence of a num- 
ber of ETs. The prospect of detecting these is limited by 
the small amount of time available per day duo to the 
Earth-Sun-ET geometry. This implies that the conven- 
tional approach to detecting asteroids by repeated obser- 
vations of a field must have cadence times in days rather 
than hours. 

Surveys of the entire ET field are impractical due to 
the observational limits imposed by the geometry, and the 
many hours it would take to conduct such a survey. This 
paper has identified the region of highest probability for 
detection, in the inclination range of 10° - 45° and helio- 
centric longitude range of 30° - 130° (L4) and 240° - 340° 
(L5). The sky area for Earth-based observers has been 
determined using numerical integrations. A strategy has 
been proposed for observing a sub-region of the ET field 
and using Earth's revolution about the Sun to progres- 
sively survey the field as it crosses this sub-region. This 
approach takes only a few minutes per session, two days 
per week, and is readily achievable by a survey telescope 
before or after the primary science mission. While this 
method requires a program of continued observations for 
one year, the total time commitment for the program is a 
few tens of hours spread throughout that year. We note 
that the only known ET, 2010 TK7, would pass through 
this region during its orbit. 

The specific observing geometry of the Gaia satellite 
and its position at Earth's L2 Lagrangian point will be 
examined in more detail in future work. Initial simula- 
tions for the detection of Trojans by Gaia show promise 
(Migiiard, pcrs. comm.). Results of detailed simulations 
will be reported with particular regard to the detection 
limits and observational mode of operation of Gaia. 

The different, and changing, geometry of the Earth- 
Mars Trojan relationship makes a detailed study more 
involved than the relatively static geometric condition of 
the ETs. As our nearest planetary neighbour with known 
Trojans, and the prediction that many more should ex- 
ist in Mars' orbit (Tabachnik & Evans 1999), a simi- 
lar modelling exercise is in progress. The peculiarities of 
the changing geometry and the implications for a search 
for additional Mars Trojans will be explored in greater 
depth. 

The search area defined for ETs may also result in 
the detection of Inner Earth Objects (lEO), asteroids 
with orbits interior to Earth's orbit about the Sun, as 
they pass across the field. To date there have been nine 
lEOs discovered however the predicted number is 36 ± 26 
with diameter > 1 km and 530 ± 240 with diameter 
> 250 m (Zavodny et al. 2008). Further work is required 



ACKNOWLEDGMENTS 

The authors would like to thank the anonymous rcforoo 
whose comments and suggestions significantly improved 
the final version of this manuscript. M. Todd thanks the 
organisers of the Gaia workshop (Pisa 2011) for providing 
a fertile environment for discussing Gaia science. D.M. 
Coward is supported by an Australian Research Council 
Future Fellowship. 



REFERENCES 

Christou A. A., Aslicr D. J., 2011, MNRAS, 414, 2965 
Connors M. ct al., 2000, BAAS, 32, 1019 
Connors M., Wiegert P.,Veillet C, 2011, Nat, 475, 481 
Drake A. J. et al., 2003, ApJ, 696, 870 
Dunbar R. S., Helin E. F., 1983, BAAS, 15, 830 
Green D. W. E., 1992, Int. Comet Q., 14, 55 
Jcdicke R., Magnier E. A., Kaiser N., Chambers K.C., 
2007, in Milani A., Valsecchi G. B., Vokrouhlick D., 
eds, Proc. lAU Symp. 236, Near Earth Objects, our 
Celestial Neighbours: Opportunity and Risk. Kluwer, 
Dordrecht, p. 341 
Jones R. L. et al., 2009, Earth Moon Planet, 105, 101 
Law N. M. et al., 2009, PASP, 121, 1395 
Mignard F. et al., 2007, Earth Moon Planet, 101, 97 
Mikkola S., Innanen K. A., 1990, AJ, 100, 290 
Morais M. H. M., MorbideUi A., 2002, Icarus, 160, 1 
Nicholson S. B., 1961, Astron. Soc. Pac. Leaflets, 8, 239 
Pilat-Lohingcr E., Dvorak R., Burge Ch., 1999, Celest. 

Mech. Dyn. Astron., 73, 117 
Scholl H., Marzari F., Tricarico P., 2005, Icarus, 175, 
397 

Tabachnik S., Evans N. W., 1999, ApJ, 517, L63 
Tabachnik S., Evans N. W., 2000a, MNRAS, 319, 63 
Tabachnik S., Evans N. W., 2000b, MNRAS, 319, 80 
Tedesco E. F., Cellino A., Zappala V., 2005, AJ, 129, 
2869 

Warner B. D., Harris A. W., Pravec P., 2009, Icarus, 

202, 134 

Wiegert P. A., Innanen K. A., Mikkola S., 1997, Nat, 
387, 686 

Whiteley R. J., Tholen D. J., 1998, Icarus, 136, 154 
Zavodny M., Jedicke R., Beshore E. C, Bernardi F., 
Larson S., 2008, Icarus, 198, 284 



© 2011 RAS, MNRAS 000, 1-5 



